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NM ABSTRACT. Simulation model of the population dynamics of N.viridula 
(Hemiptera: Pentatomidae). In the north of the Buenos Aires province, the 
green stink bug Nezara viridula develops 3 almost discrete generations along 
an activity period which lasts approximately 28 weeks. An intensive study du- 
ring 9 generations allowed the construction of life tables and a key factor 
analysis was performed. N. viridula suffers the parasitism of eggs by the sce- 
lionid Trissolcus basalis and nymphal and adult predation and parasitism by 
the tachinid Trichopoda giacomellií. The tachinid reduces in a density depen- 
dent manner the potencial fecundity because of the reduction of adult longe- 
vity, without affecting egg fertility. In this work, a preliminar projection matrix 
model of the population dynamics of N. viridula is presented. The time inter- 
val is the week, and the following developmental stages (with the correspon- 
ding developmental time “rounded” to week units) were considered: eggs (1 
week); nymphs l-I (2 weeks); nymps IV (1 week); nymphs V (2 weeks) and 
adults (17 weeks, maximum longevity). The survivorship of the eggs and 
nymphs were considerd as constat during a given generation but changes bet- 
ween generations according to the key factor analysis. The numerical respon- 
se of T. giacomellii introduces a temporal delay of 5 weeks in the expression 
of the mortality of adult hosts. The model mimics satisfactorily the density of 
eggs, nymphs and adults observed, and the effect of T. giacomellii is capable 
of regulate the system even if the time lags and the exponent that represent the 
intensity of adult parasitism were varied 40 %. 


INTRODUCCIÓN 


La finalidad del desarrollo de modelos teóricos 
de las interacciones entre presas y predadores así 





como entre huéspedes y parasitoides ha sido el 
de acrecentar nuestro conocimiento sobre estos 
sistemas, para, entre otras posibles metas, colabo- 
rar en el manejo de plagas en distintos agroeco- 
sistemas. Si bien estos objetivos no han sido com- 
pletamente alcanzados aún, diversos enfoques y 
desarrollos teóricos permiten ser optimistas. La 
sugerencia hacia el desarrollo de modelos con 
suficiente detalle del sistema en estudio, que si 
bien pierden la generalidad de los modelos teóri- 
cos clásicos ganan en realismo, ha sido uno de 
aquellos desarrollos (Godfray & Waage, 1991; 
Murdoch & Briggs, 1996). En particular, los mo- 
delos matriciales discretos con consideración de 
la estructura etaria de la población, debido a su 
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flexibilidad han sido extensamente usados tanto 
en el estudio de la biología y dinámica poblacio- 
nal, así como en las interacción de competencia 
y predación de diversas especies de plantas y ani- 
males. Estos modelos, proyectan temporalmente 
el estado de una población estructurada, inte- 
grando la información sobre supervivencia, fe- 
cundidad, desarrollo y crecimiento específicos 
(Groenendael et al., 1988). 

Nezara viridula es un hemíptero polifitófago 
de distribución cosmopolita, y en la presentación 
(“Análisis del parasitismo en Nezara viridula por 
los parasitoides Trissolcus basalis y Trichopoda 
giacomellii en términos de la reducción de pará- 
metros demográficos”) se detallan las característi- 
cas del ciclo de vida así como los resultados del 
análisis de tablas de vida en laboratorio y a cam- 
po. Se describen asimismo las características den- 
sodependientes del parasitismo por las especies 
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mencionadas, en particular, la capacidad de re- 
gulación de T. giacomellii. 

En este simposio se presenta un modelo de si- 
mulación preliminar de la dinámica poblacional 
de N. viridula, considerado la estructura básica 
para desarrollos posteriores, dejando para otras 
contribuciones el análisis del modelo. La finali- 
dad es ayudar tanto a interpretar la dinámica del 
sistema huésped-parasitoide como a la toma de 
decisiones para el manejo de N. viridula. 

El modelo de la interacción. El modelo de la 
dinámica poblacional de N. viridula es un “Pro- 
jection matrix model” (Groenendael et al., 1988), 
consistente en una matriz de Leslie (Leslie, 1954) 
modificada, ya que los valores de las distribucio- 
nes de fecundidad y supervivencia específicos no 
son constantes, sino que varían en las sucesivas 
generaciones, o bien son funciones que represen- 
tan una densodependencia retrasada. El mismo 
fue escrito en lenguaje FORTRAN, y esquemáti- 
camente puede representarse como en la tabla 1. 

De acuerdo con la periodicidad semanal de los 
muestreos desarrollados en el estudio poblacional 
en el campo, el modelo “corre” también a interva- 
los (t) semanales, hasta completar las aproximada- 
mente 28 semanas que dura un período de activi- 
dad, en el que se suceden tres generaciones. En el 
modelo se considera explícitamente la estructura 
de edades (j), abarcando individuos de 1 hasta de 
24 semanas de vida. La correspondencia entre la 
edad del individuo y su estado de desarrollo es la 
siguiente: huevo (1 semana: ¡=1); ninfas |-111 (2 se- 
manas: ¡=2, j=3); ninfas IV (1 semana: j=4); ninfas 
V (2 semanas: j=5, j=6) y el adulto (18 semanas: 
j=7,....., j=24). Esta correspondencia se efectuó al 


“redondear” en semanas, la duración de los esta- 
dos de desarrollo estimada en días y en condicio- 
nes de laboratorio a 22 °C (la temperatura prome- 
dio del período de actividad) (Liljesthróm, 1983). 

Las distribuciones de fecundidad específica, 
m(¡,h), y supervivencia específicas, s(¡,h), ocupan 
la primera fila y la subdiagonal principal, respec- 
tivamente, y varían según la generación (h) de que 
se trate: la Tra. generación dura ocho semanas 
(t=1 a t=8), la 2da. y la 3ra., 10 semanas cada una 
(t=9 a t=18) y (t=19 a t=28), respectivamente. 

La distribución de fecundidad de la primera ge- 
neración se calculó en condiciones de laboratorio, 
correspondiendo en consecuencia la notación del 
modelo a m(¡,1). Además, el período pre-oviposi- 
cional dura alrededor de dos semanas, por lo cual 
el valor de fecundidad para la segunda semana es 
nulo: m(8, 1)=0. Las siguientes generaciones de N. 
viridula no varían respecto al número promedio de 
posturas que una hembra efectuaría durante su vi- 
da, aunque sí exhiben un aumento en al tamaño 
promedio de las posturas: 66,87 huevos la primera 
generación; 79,27 huevos la segunda y 79,92 hue- 
vos la tercera. En consecuencia, para representar 
las distribuciones de fecundidad correspondientes 
a la segunda y tercera generaciones, se multiplicó 
la distribución de fecundidad específica estimada 
más arriba por un factor R(h), que es el cociente 
entre el tamaño promedio de las posturas de las su- 
cesivas generaciones y 66,87. 

La supervivencia de huevos y ninfas de cada 
generación se obtuvo a partir del análisis de la 
mortalidad durante nueve generaciones sucesivas 
(Liljesthróm & Bernstein, 1990) y la supervivencia 
de adultos se considera densodependiente retra- 





Tabla I: Estado de desarrollo, duración y supervivencia de N. viridula. 


Estado de desarrollo 


en la generación h 


Duración Supervivencia s(j, h) 


(edad y generación h) 


as p» huevos 1 semana 

: A > s(1,h) 

i Ninfas I-II 2 semanas 

| E R E æ  s(2,h) y s(3,h) 
| Ninfas IV 1 semana 

0 E >æ s(4h) 

i NInfas V 2 semanas 
a æ s(5,h) y s(6,h) 
; Adulto prereproductor l a 2 semanas: f(t) 

E e E æ s(7,h); s(8,h)p 
i Adulto reproductor 18 semanas 

A O =--> s(9,h),....s(24,h) 
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sada (Liljesthróm, 1993), y en ambos casos fueron 
expresadas en términos del factor k = log densi- 
dad inicial/ log densidad final (Varley € Grad- 
well, 1963). Simbólicamente: n (t+1) = An (0; 
donde: n(0) = (nl, n2, n3,........ ,n24)(t) y n (t+1) = 
A E ,n24)(t+1} son vectores columna 
que representan los individuos de edad j=1 hasta 
j=24 de N. viridula en las semanas (t ) y (t+1); A 
es una matriz cuadrada de 24 x 24 (Tabla II): don- 
de, m(¡,h) representa la distribución de fecundi- 
dad específica calculada en el laboratorio y R(1) 
= 1; R(2) = 1,185 y RGB) = 1,195. S(1,h) = 10 exp 
[- k(2),h] representa (respecto de la generación h- 
ésima) la supervivencia de huevos al parasitismo 
por T. basalis fundamentalmente. 


Tabla II: Matriz cuadrada de 24x24 (explicación en el 
texto) 


0 0 0 0 0 0 0 R(h)mGZh)..R(h) m(4h) 
sah) 0 0 0 0 0O 0 Dad 0 
0SQH)0 0 0 0 0 O. T 
0 osh) o 0 0 0 ata 0 
0 0 0S(4h)0 0 0 0... ad 
Pa 0 0 05S56ho0 0 OEE 0 
0 0 0 0 O0S(6h)O0 E TS 0 
TERRE E IA 30 
A E E T O 0 
0 0 0 0 0 0 0 S(23L() 0 
La mortalidad de las ninfas l-IH (predación y 


efectos de factores climáticos) se concentra en la 
primera semana (siendo nula en la segunda), por lo 
tanto, las respectivas supervivencias se represen- 
tan: SQ,h) =10 exp I- k(3),h]; y S(3,h)=1. La super- 
vivencia de las ninfas IV y ninfas V (predación y pa- 
rasitismo por T. giacomelli¡) se calcularon de mane- 
ra similar a las anteriores, siendo sus respectivos 
valores: S(4,h) = exp [k(4)h]; y S(5,h) = exp [-k(5),h] 
(como antes, dado que en las ninfas V toda la mor- 
talidad se concentró en la primer semana que dura 
el estadio, la supervivencia correspondiente a la se- 
gunda semana es la unidad: S(6,h) = 1). 

Los valores de k(2), k(3), k(4) y k(5) represen- 
tan la intensidad de mortalidad en términos del 
factor k, y se encuentran en la Tabla 1 de Lil- 
jestrthóm & Bernstein (1990). 

La supervivencia de los adultos debido al pa- 
rasitismo por T. giacomellii : S(j) (toma los valores 
estimados en condiciones de laboratorio (prome- 


67 


dio de 5 tablas de vida en el apéndice de Liljesth- 
róm, 1983) desde j=7 a j=24, y la función L(t) = 1 
para todo j < LM+S; y L(t) =0 para todo j > LM+S; 
siendo: LM = 28,45 [[ E n (j, t-0) | exp (-0,247)), y 
S las semanas que dura el desarrollo de huevo a 
adulto de N.viridula. 

La mortalidad de los adultos (debida funda- 
mentalmente al parasitismo por T. giacomellii, 
que no exhibe preferencia respecto de la edad de 
los huéspedes adultos, se manifiesta por la reduc- 
ción de la longevidad potencial (limitando así la 
fecundidad), pero sin afectar significativamente la 
fertilidad. El efecto es claramente densodepen- 
diente con retraso temporal, y está expresado en 
términos del factor k por la siguiente ecuación de 
regresión lineal: log (LP/LM) = 0,43 + 0,247 log 
A; donde LP = 76,54 días representa la longevi- 
dad potencial (Liljesthróm, 1983), LM representa 
la longevidad generacional media en condiciones 
de campo, y A representa la densidad generacio- 
nal de adultos de N. viridula (Liljesthróm, 1993). 
Luego, despejando y expresando la densidad ge- 
neracional A = Y n(¡,h), queda :LM=28,45 {[ Y 
n(,h), exp €0,247)) . Sin embargo, debido al cla- 
ro retraso temporal (t) en el efecto del parasitoi- 
de, se tomó la densidad de adultos correspon- 
diente a t=5 Semanas antes, tiempo que represen- 
ta la duración aproximada del desarrollo de hue- 
vo a adulto de T. gracomellii. 


RESULTADOS Y DISCUSIÓN 


En las figuras 1 y 2 se muestran las densidades 


observadas y calculadas, respectivamente, de nin- 


fas IV, ninfas V y adultos de N. viridula correspon- 
dientes al período de actividad 1981-1982. En la fi- 
gura 3 se muestran, para el mismo período, las den- 
sidades observadas y calculadas de huevos; y en la 
figura 4, las densidades observadas de adultos (co- 
rrespondientes a los períodos de actividad 1981- 
1982; 1982-1983 y 1983-1984) y calculadas. 

Se observan tres generaciones casi discretas, 
siendo la segunda generación de adultos de cada 
período la más densa. Estos picos, sin embargo, 
muestran una considerable variación, en función 
de la densidad inicial de adultos postinvernantes 
y las mortalidades de los estados preimaginales, 
en particular la que afecta a las ninfas 1-11, que en 
resultados anteriores fue identificada como el fac- 
tor clave en la dinámica poblacional de N. viridu- 
la. Es de destacar, respecto de los adultos de esta 
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segunda generación, que son los que comenza- 
rán a invadir cultivos de soja cercanos. 


Individuos/m? 
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Fig. 1. Variaciones numéricas observadas de ninfas IV, 
ninfas V y adultos de N.viridula, correspondiente al 
período de actividad 1981-1982. 
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Fig. 2. Variaciones numéricas calculadas de ninfas IV, 
ninfas V y adultos de N.viridula, correspondiente al pe- 
ríodo de actividad 1981-1982, con los siguientes valores 
iniciales: adultos prereproductores, n(7,1)=2,74; la su- 
pervivencia de huevos, ninfas I-HI, ninfas IV y ninfas V, 
en la primera generación: s(1,1)=0,62; s(2,1)=0,56; 
s(3,D=1; s(4,1)=0,8; s(5,1)=1; s(6,1)=0,94. En la segunda 
generación: s(1,2)=0,19; s(2,2)=0,06; s(3,2)=1; s(4,2)=0,95; 
s(5,2)=1; s(6,2)=0,78; y en la tercera generación: 
s(1,3)=0,43; s(2,3)=0,15; s(3,3)=1; s(4,3)=0,88; s(5,3)=1; y 
s(6,3)=0,83. El retraso temporal fue: t = 5 semanas. 
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Fig. 3. Variaciones numéricas observadas y calculadas 
de huevos de N.viridula, para el período de actividad 
1981-1982 (los valores iniciales para los calculados, son 
los mismos que en la Fig. 2). 
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Fig. 4. Variaciones numéricas observadas y calculadas 
de adultos de N.viridula, correspondiente a los perío- 
dos de actividad 1981-82; 1982-83 y 1983-84. Los valo- 
res calculados parten de la siguientes condiciones ini- 
ciales: adultos preimaginales, n(7,1)= 2,74; supervi- 
vencia de huevos, ninfas l-JII, ninfas IV y ninfas V, pa- 
ra las tres generaciones: s(1,)=0,43; s(2,))=0,16; s(3,j)=1; 
s(4,)0,88; s(53)=1; y s(6,)= 0,83. El retraso temporal 
fue: t = 5 semanas. 


Los valores calculados describen aceptable- 
mente las variaciones numéricas, aunque se debe 
recordar que la información obtenida de estos 
muestreos en el campo, fue usada para generar el 
modelo. Las variaciones numéricas de adultos 
observadas en otros cuatro períodos de actividad 
(de los que se desconocen estimaciones de super- 
vivencia preimaginal) exhibieron un patrón simi- 
lar, tal como se muestra en la Fig. 4 en Liljesthróm 
& Bernstein (1990) y en las Fig. Ta y 1b en Lil- 
jesthróm £ Cameán (1992). Se diferencian, como 
antes, en los picos alcanzados por la segunda ge- 
neración. 

Como un dato preliminar respecto del com- 
portamiento del modelo, se muestra en la figura 5 
las variaciones numéricas de adultos de N. viri- 
dula, partiendo de idénticas condiciones iniciales 
a los de la figura 4. La “corrida” representa apro- 
ximadamente ocho períodos de actividad conse- 
cutivos, sin el efecto del “corte” invernal que im- 
plicaría: 1) un reordenamiento de la estructura de 
edades al comienzo de cada período (sólo existen 
adultos prereproductores) y 2) una marcada re- 
ducción del parasitismo por T. giacomellii al co- 
mienzo de cada período por una baja densidad 
inicial de parasitoides adultos. Se observa que la 
densidad final fluctúa en un rango estrecho. La 
persistencia y relativamente alta estabilidad del 
sistema se debe, tal como se puntualiza en resul- 
tados anteriores, al parasitismo densodependien- 
te por T. giacomellii, y el retraso temporal intro- 
ducido en el modelo, provoca las oscilaciones en 
la densidad. 
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Como es conocido, al aumentar la magnitud 
del retraso temporal, se aumenta la amplitud de 
las variaciones de densidad: en las figuras 6 y 7 
se muestran, bajo la mismas condiciones iniciales 
que en la figura 4, las “corridas” en las que se in- 
trodujo una modificación en la magnitud del re- 
traso temporal, En la figura 6, el mismo se redujo 
en 40% (de 5 días a sólo 3 días), y en la figura 7 
se aumentó en 40% (de 5 días a 7 días), variando 
ampliamente, como se observa, el rango de osci- 
lación de la densidad. 

Respecto de la densodependencia introduci- 
da en el modelo, cuando la misma desaparece 
(el coeficiente m = 0 en la ecuación: LM = 28,45 
(E n (j, t-t) | exp (=m)), y la mortalidad resultan- 
te es una constante), la población desaparece, ya 
que el conjunto de las mortalidades supera a la 
fecundidad. En las figuras 8 y 9 se muestran, res- 
pectivamente resultados obtenidos al aumentar 
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Fig. 5. Variaciones numéricas calculadas de adultos de 
N.viridula, a partir de las mismas condiciones iniciales 
que en la Fig. 4. 
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Fig. 6. Idem Fig. 4, excepto que el retraso temporal es 
de 3 días. 


el valor de m en 40% (de m=0,247 a m=0,3458) 
y a disminuirlo en 40% (de m=0,247 a m=0,1482). 
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Fig. 7. Idem Fig.4, excepto que el retraso temporal es 
de 7 días. 
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Fig. 8. Idem Fig. 4, excepto que el coeficiente m varía 
de m=0,247 a m=0,3458. 
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Fig. 9. Idem Fig. 4, excepto que el coeficiente m varía 
de m=0,247 a m=0,1482. 
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